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Abstract 

We calculate total and differential production rates of D* mesons in two-photon colli- 
sions at LEP2. We include full next-to-leading order QCD corrections, and perform an 
extensive study of the sensitivity of our predictions to variations of the renormalization 
scale, charm mass, photonic parton distribution set and fragmentation function. The 
results are compared with recent data from LEP2. 
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1. INTRODUCTION 



Charm quark production in two-photon collisions at high-energy e + e~ colliders provides 
new possibilities to study the dynamics of heavy quark production; it complements the 
extensive analyses that have been carried out at fixed-target experiments and at other 
high-energy colliders jy]. In two-photon collisions^] each of the photons can behave as either 
a point-like or a hadronic particle. Consequently, one distinguishes in such collisions direct- 
(both photons are point-like), single resolved- (one photon is point-like, the other hadron- 
like), and double resolved (both photons are hadron-like) production channels. The resolved 
channels require the use of parton densities in the photon, whereas the production via the 
direct channel is free of such phenomenological inputs. In general, the different channels 
mix when including higher orders in perturbation theory, and thus the distinction between 
the direct and resolved contributions becomes non-physical and scheme-dependent. 

The mass of the heavy quark, m Q Aqcd, acts as a collinear cutoff and sets the 
hard scale for the perturbative calculation at small heavy-quark transverse momentum p T . 
It is thus possible to define an all-order infrared-safe cross section for open heavy flavour 
production, even at p T = 0. The heavy quark mass also ensures that the separation into 
direct and resolved production channels is unambiguous at next-to- leading order (NLO). In 
the case of charm production, however, the heavy quark mass is not very large with respect 
to Aqcd, and one therefore expects large radiative corrections in perturbative QCD. This 
is indeed borne out by NLO analyses of charm production in hadron-hadron and photon- 
hadron collisions, where the leading-order result sometimes only accounts for less than 50% 
of the full next-to-leading order rate JIJ. As will be discussed at length in what follows, the 
situation in photon-photon collisions is quite different, and the predictions are under better 
theoretical control, in spite of the presence of resolved channels. 

The OPAL collaboration has presented 0] new data for D* production in two-photon 
collisions, at (mostly) yfs e + e - = 189 GeV. Besides the total cross section cr^ , OPAL has 
measured the differential rate with respect to the D* transverse momentum, da®* /dp®* , 
and pseudorapidity da®^ /dr] D . Recently, also L3 has presented results for D* production 
in two-photon collisions [||. L3, however, does not apply an antitag condition for the 
scattered electrons. It is thus not obvious how the corresponding data can be compared 
with a theoretical calculation where the cross section is defined in the Weizsacker- Williams 
approximation in the standard way. 

The new LEP2 data motivated us to perform an extensive study of D* production in two- 
photon collisions, including all NLO QCD corrections. We use fully differential NLO Monte 
Carlo programs for all the production mechanisms, constructed in refs. 0-^]. Earlier studies 
either focused on an individual component of the production Q , or had mostly HERA data 
in mind fill, or did not include fragmentation functions 0. Moreover, a relatively large 
number of NLO-evolved photonic parton distributions [§-|Tl|l, and new studies on charm 
fragmentation functions [12-14] are available, some of them incorporating the most recent 
theoretical and experimental findings. We shall examine in detail the sensitivity of the NLO 



1 In this paper, we will only consider the case in which the incoming photons are on-shell. 
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calculations of D* meson cross sections with respect to variations of the renormalization 
scale, charm mass, photonic parton distribution sets, fragmentation function, and others. 

The paper is organized as follows. In section 2 we review very briefly the NLO Monte 
Carlo programs that we used in the rest of the paper. The results of our study, and 
the comparison of our calculations with the OPAL data, are presented in section 3. We 
summarize our work in section 4. The computation of photon-photon cross sections and the 
issue of scheme dependence of the direct and resolved production mechanisms is discussed 
in the appendix. 



2. NLO MONTE CARLO PROGRAMS 

We consider the process 

e + + e~ — ► e + + e~ + Q + Q, (1) 

where the final state positron and electron are scattered almost collinearly to the beam line. 
The cross section is dominated by the scattering of two on-shell photons, which produce 
the heavy quarks and can be written as a convolution of two Weizsacker- Williams spectra 



[15] with the cross section for the process 77 — ► QQ. As mentioned in the introduction, 
besides the direct contribution, also resolved production channels have to be taken into 
account. The cross section for the resolved contributions factorizes into a partonic hard- 
scattering cross section, convoluted with photonic light-quark and gluon densities. Notice 
that these densities grow as a em /a s at large scales, owing to the inhomogeneous term in the 
relevant Altarelli-Parisi equations; therefore, the cross sections for the resolved processes 
are formally of the same perturbative order as the cross section for the direct process. More 
details on the computation of photon-photon rates in perturbative QCD will be given in 
the appendix. 

To calculate the direct and resolved contributions including the complete NLO correc- 
tions, we use Monte Carlo programs that are fully exclusive in all final-state particles. We 
want to emphasize that our programs are different from the usual Monte-Carlo parton- 
shower codes, since they result from fixed-order QCD computations, where no collinear 
approximation (and no shower evolution) has been performed. For the direct channel, the 
NLO program has been constructed in Q , while for the single-resolved and double-resolved 
channels we use the NLO programs constructed in refs. j|,|5|]. To check our results we 
employed in addition the single-particle inclusive programs of refs. |7|,[l6.17]. In our codes 



all final-state kinematical quantities are available on an event-by-event basis, and it is thus 
possible to calculate any infrared-safe quantity to NLO, if the quantity requires two or three 
partons for its definition, or to LO if it requires three partons, with any final-state kine- 
matical cuts matching those implemented by experiments. This allows us to compare NLO 
predictions with measurements in the experimentally visible region. As will be discussed 
below, cross section extrapolation beyond the accepted region introduces large theoretical 
uncertainties, at least for the case at hand, and we find the visible cross section to be as 
important as the total cross section for the comparison of data with theory. The wide avail- 
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ability of flexible higher order Monte Carlo programs in which acceptance cuts can be easily 
built in makes the calculation of many types of visible cross section practical, thereby pro- 
viding the possibility of a more detailed comparison between theoretical and experimental 
results at high energy colliders. 



3. RESULTS 

We will first define a default set of input parameters and then study parameter varia- 
tions to examine the sensitivity of our theoretical predictions. The collision energy is fixed 
at ^/s e + e - = 189 GeV. We set the charm quark mass, defined in the on-shell renormaliza- 
tion scheme, to m c = 1.5 GeV. The two-loop expression for the QCD coupling is used, with 

(5) 

A— = 220 MeV, as favoured by recent global QCD fits. The default choice for the renor- 

MS ' 

malization scale is p? R = m% = m^+p T (c) 2 , and for the factorization scale we use p? F = Am^,. 
Our default choice for the photonic parton densities entering the resolved channels is the re- 
cent NLO-evolved set GRS-HO 0. We checked that, by adopting the GRS-HO set defined 
in the MS or DIS 7 scheme, and changing the short-distance cross sections accordingly, 
only numerically negligible differences arise. Finally, for the Weizsacker-Williams function 
we adopt the form proposed in ref. [O], designed to describe an antitag condition on the 
scattered electron. Here, the antitag angle is # max = 0.033, as determined by the OPAL 
experimental setup. The Peterson et al. parametrization |2(J is used as the charm-to-Z)* 
fragmentation function, normalized according to Jq 1 dzD(z,e) = 1. Our default value for 
the non-perturbative parameter entering the fragmentation function will be e = 0.035 (see 
for example ref. flUf ). In order to obtain the D* kinematical variables from those of the 
fragmenting charm quark, we rescale the three space components of the parton momentum 
with the momentum fraction z, and compute the energy according to the charm quark mass 
shell condition (the choice of mass to fix the D* energy component is arbitrary, because the 
factorization theorem only holds for large p T , where mass effects are negligible. We have 
checked that our numerical results in the experimentally visible region p T (D*) ^> 2 GeV are 
not sensitive to the specific choice). 

We have to comment on the fact that the Aqcd value we use is not the one associated 
with the GRS-HO photon set. The parton densities are correlated with Aqcd Q> so that, in 
principle, Aqcd cannot be chosen independently of the PDF set. However, in practice, it is 



known that the correlation between Aqcd and the densities of the photon is very mild [21], 
and our choice can be considered safe, also given the fact that other theoretical uncertainties 
are dominant in charm physics. Since our aim is to compare theory with experimental data 
and because the normalization of the resolved contributions is very sensitive to the choice 
of Aqcd, we deem it appropriate to use an up-to-date value for Aqcd- As a final comment 
on our default choices, we remark that we prefer to use \x F = 2m T rather than p, F = m T as 
the factorization scale: at very low transverse momenta, fj, F = m T would probe the parton 
densities at a scale where the validity of the evolution equations is not on firm grounds. As 



2 The correlation is due to the Altarelli-Parisi evolution equation, and to the fact that the largest sensi- 
tivity is to asg(x), where g(x) is the gluon density, and not to g(x). 
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is customary in charm physics, we will not consider the variation of our results with respect 
to fi F , for the same reason. The scale uncertainty that we will quote therefore only partially 
accounts for the full theoretical error. 

3.1. TOTAL CROSS SECTIONS 

The preliminary value obtained by OPAL for the total cross section for D* production in 
77 collisions at ^Js e + e - = 189 GeV, in the range 2 < p T (D*) < 12 GeV and with the cut 
\t](D*)\ < 1.5 is 

a(e + e~ - e+ e~ D* X) [2<Pt{d , )<12 GeV ;|„(D*)|<i.5] = 29.9 ± 4.2 pb (OPAL). (2) 

Taking into account the probability for a charm quark to fragment into a D* meson (we 
use here f(c — ► D* + ) = 0.233 ± 0.010, the same value as adopted in ref. 0), and the fact 
that OPAL counts both the heavy quark and the heavy antiquark in their cross-section 
determination, we find a theoretical value of 

a(e + e~ e + e~ D* X) [2<Pt (d*)<\2 GeV;KD*)[<i.5] = 17 - 3 ^9 pb ( NL0 Q CD )- ( 3 ) 

The theoretical error reflects the dependence on the charm quark mass, chosen in the range 
1.2 < m c < 1.8 GeV, and the renormalization scale uncertainty, m T /2 < fi R < 2m T . The 
renormalization scale has been varied independently for the direct and resolved contribu- 
tions. We observe that the central NLO prediction underestimates the experimental data, 
which seems to suggest that the input parameters should be taken such as to yield the largest 
possible cross section. This issue will be discussed in more detail in the next subsection. 

To provide more detailed information on the NLO QCD cross section prediction and on 
the theoretical uncertainties, we have collected the results for the total and the visible cross 
section in table [I], including the variation of different input parameters. We consider the 
effect of changing the renormalization scale and the charm quark mass, and of adopting the 
photon density sets GRV-HO || and AFG-HO |J. We also show the sensitivity to Aqcd 
by choosing A^ = 200 MeV, which is equal to the value associated with GRV-HO and is 
smaller than that associated with GRS-HO. The first two columns of the table present total 
cross sections in the full phase space (i.e. no final-state kinematical cuts are applied); they 
differ in the choice of the reference scale /j,q, which in the first column is set to a fixed value 
(the charm mass), while in the second column it depends on the final-state kinematics (the 
scale is set equal to the transverse mass). Since the total cross section is dominated by small 
transverse momenta, p T ~ m c , the difference between choosing a fixed or a p T -dependent 
scale is less than about 10% for the sum of direct and resolved contributions. Finally, the 
third column displays the predictions for the visible range, defined by the OPAL cuts. It 
is obvious from table [l] that the visible cross section is much more stable under scale and 
mass variations than the fully extrapolated one and thus provides an important and in 
certain aspects better quantity for comparison of data with theory. The dependence of the 
theoretical uncertainties upon the kinematical cuts will also be studied in the following. 
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Mo = m T 



°DD* (P b ) 




Mo = ra c 


Mo = rn T 


2 < p T < 12 GeV 
|?y| <, i.o, e — u.uoo 


Total 


Default 


871.1 


827.8 


37.22 




Mb = 2/xo 


742.2 


727.1 


35.46 




H R = no/2 


1204 


1035 


40.45 




m c = 1.2 GeV 


1712 


1591 


41.70 




m c = 1.8 GeV 


504.8 


485.6 


32.7 




A( 4 ) = 200 MeV 


751.2 


733.0 


35.35 




GRV-HO 


1053 


984.6 


42.14 




AFG-HO 


906.4 


856.2 


39.74 


Direct 


Default 


590.5 


577.5 


26.20 




Ma = 2/U 


528.9 


522.8 


26.60 




Ms = Mo/2 


772.7 


723.0 


25.57 




m c = 1.2 GeV 


1098 


1066 


27.79 




m c = 1.8 GeV 


357.3 


350.9 


24.10 




A^ 4 ) = 200 MeV 


546.0 


538.4 


26.46 


Resolved 


Default 


280.6 


250.3 


11.02 




Ma = 2/x 


213.3 


204.3 


8.860 




Ms = Mo/2 


430.8 


311.5 


14.88 




m c = 1.2 GeV 


613.9 


524.5 


13.91 




m c = 1.8 GeV 


147.5 


134.7 


8.659 




A( 4 ) = 200 MeV 


205.2 


194.6 


8.894 




GRV-HO 


462.6 


407.1 


15.94 




AFG-HO 


315.9 


278.7 


13.54 



Table 1: Results for the total cross section broken down to the direct and the (single- 
plus double- ) resolved photon contributions. The results in the first two columns are 
independent of the value of e. The factorization scale is fi F = 2^ - The renormal- 
ization scale is fj, R = Hq, unless otherwise indicated. Note that the fragmentation 
function is normalized according to J Q dzD(z, e) = 1 as described in the text. 
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3.2. DIFFERENTIAL DISTRIBUTIONS 



The study of differential distributions constitutes a more comprehensive test of the theory 
than that of total cross sections alone. In this section we concentrate on the single-inclusive 
distributions in transverse momentum p T and pseudorapidity rj, for which experimental data 
are now available. We emphasize that with our Monte Carlo programs also more exclusive 
observables can be calculated at NLO, so that additional comparisons can be made, should 
experimental data become available in the future. In the following, we have defined the Born 
cross section for the resolved channels as the lowest-order partonic cross sections computed 
with two-loop a s and convoluted with NLO-evolved photonic parton density sets. This 
choice implies that the difference between the Born and NLO results is essentially due to 
the short-distance parton dynamics, as is the case of the direct part. 

In fig. [l] we show our predictions for the p T and rj distributions using the default set 
of parameters. The cuts \rj(D*)\ < 1.5 and 2 < p T (D*) < 15 GeV have been applied, 
respectively. The total cross sections are decomposed into the direct, single-resolved and 
double-resolved components. We stress again that these quantities are not physical in 
general, although at this order the direct component is well defined as such. The relative size 
of the single- and double-resolved components is factorization-scheme- and scale-dependent; 
however, for both the MS and DIS 7 schemes, and regardless of the specific choice for ji F , 
the double-resolved component contributes only very marginally to the charm cross section 
at LEP2 energies. The direct channel dominates the cross section in the visible range, 
even more so for large p T . As for the rj distribution, one observes that the direct channel 
dominates at central values, but that the single-resolved contribution extends further out in 
pseudorapidity. Therefore, when only large values of rj are considered, there is an enhanced 
sensitivity to resolved contributions, making it in principle a good observable for learning 
about the partonic densities 0. However, in practice it is very hard for the LEP experiments 
to properly identify charmed mesons in such forward regions. Therefore, we shall not pursue 
the study of such observables here. 

Let us briefly comment on the effects of resumming large logarithms log(p T /m c ), which 
arise beyond LO from the collinear emission of gluons by heavy quarks at large p T , or from 
almost collinear branching of gluons into charm-quark pairs. At p T 3> m c these terms might 
spoil the convergence of the perturbation series and lead to large scale dependences of the 
NLO result. The large logarithms can be resummed by neglecting in the matrix elements all 
the power-suppressed mass terms (which behave like m c /p T to some power; this is the reason 
why this approach is often - not very rigorously - denoted as "massless"), and absorbing 
the mass singularities, occurring as powers of log(p T /m c ), into parton distribution functions 
and parton-to-heavy-quark fragmentation functions. The fragmentation functions can be 
computed in perturbation theory [^] at a scale of the order of the heavy quark mass, and 
subsequently evolved to the appropriate scale (which is of the order of p T ) by using the 
Altarelli-Parisi equations; in this evolution, large logarithms are properly resummed. A 
study of such a resummation in photon-photon collisions has been presented in ref. [23]. 

3 Similar observables for inclusive jet final states at the Tevatron are indeed used for such purposes. 
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Figure 1: NLO predictions for the D* transverse momentum and pseudorapidity 
distributions, with default parameters (solid line). The direct, single-resolved (1-res) 
and double-resolved (2-res) results are also shown. 



We would like to emphasize that the resummed calculation is not reliable unless p T S> m c , 
since all non-singular mass terms are neglected. We have verified numerically that omitting 
those mass terms from the NLO cross section overestimates the full result by no less than 



« 100% at m c /p T « 1 and still by « 25% at m c /p T 1/4. The observation of ref. [23], 
that resummed and fixed-order calculation agree down to relatively small values of p T , must 
therefore be considered accidental. Similar findings have been made in an analysis of the 



p T spectrum in heavy-flavour hadroproduction 24]. We conclude that a comparison of the 
photon-photon D* data and the resummed approach is not meaningful unless Pt(c)<; 6- 
7 GeV. 

The relative contribution of the direct, single- and double-resolved components to the 
differential cross section is displayed in fig. | as a function of p T and r\. The curves are 
shown for both Born and NLO cross sections. In most of the p T range the relative fractions 
at the Born and NLO level are rather similar, even more so for large p T , since the value 
of a s (m T ) is decreasing with increasing transverse momentum. On the other hand, the 
relative fraction changes sizeably over the whole range of r\ when radiative corrections are 
included, because the visible cross section is dominated by the region of low p T <; 2 GeV. In 
general, the radiative corrections tend to increase the importance of the resolved processes. 
Although the relative contribution of the resolved cross section never exceeds 30% (with the 
default set of parameters), NNLO corrections may increase the importance of the resolved 
channels and thereby enhance the total D* cross section. The NLO single-resolved fraction 
falls below the LO value at low p T <, 1 GeV, causing the NLO direct fraction to exceed the 
LO one. This is due to the jq subprocess in the single-resolved component, which is only 
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e e — > e e D X: 
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1-res (NLO) 
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2- res (NLO) 
2-res (LO) 



2 4 6 
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P T [GeV] 
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0.8 - 



e + e — > e + e~ D X: 

drj/dr|(direct/l-res/2-res) / do/dr|(total) 

Vs = 189 GeV; 2 GeV < p T (D*) < 15 GeV 



0.6 
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0.2 



direct (NLO) 
direct (LO) 

1-res (NLO) 
1-res (LO) 



0.5 



1.5 



2-res (NLO) 
2-res (LO) 



2.5 3 3.5 
pseudorapidity r\ 



Figure 2: Relative contributions of the direct, single-resolved and double-resolved 
processes to the total cross section, as a function of transverse momentum and pseu- 
dorapidity. Both the NLO (histograms) and Born (symbols) results are presented. 
Default parameters as specified in the text. 



present from NLO onwards and which gives a negative contribution at low p T . The same 
behaviour turns out to occur for the double-resolved component, where the ^-initiated 
subprocess decreases the NLO cross section at very small p T . 

In fig. [H we show the ratio of the NLO cross section over the Born cross section (K- 
factor), as a function of p T and r\. The if-factor is decreasing with increasing transverse 
momentum, a consequence of the decreasing a s (m T ) and the increasing relative size of the 
direct contribution. In the visible range 2 < p T < 15 GeV the ET-factor is moderate, ranging 
approximately from 1.2 to 0.75. This has to be compared with the case of photon-hadron 
and hadron-hadron collisions, where much larger if- factors (up to K ~ 2) are observed and 
where NNLO corrections are expected to be sizeable. The dependence of the if-factor on 
the pseudorapidity r/, also shown in ^, is very mild. 

Let us now turn to the sensitivity of the transverse momentum distribution to variations 
of the input parameters. Figure. |] shows the p T spectrum of the D* for different choices for 
the renormalization scale, the charm quark mass, the NLO photonic parton density set, and 
the shape of the Peterson et al. fragmentation function as controlled by the parameter e. 
The scale, mass and parton density dependence is sizeable at small p T <; 2 GeV, the biggest 
uncertainty coming from the variation of the charm quark mass. The strong renormalization 
scale and charm mass dependence is mainly induced by the large variation in the value of 
a s [m T ). For the lowest charm quark mass value considered, and at small p T , the hard scale 
of the process, m T ~ 1.5 GeV, may indeed be too small to allow for a reliable perturbative 
analysis. However, the situation is much improved in the region probed by experiment, 



8 



2 



2 




pseudorapidity r\ 



Figure 3: Ratio of NLO over Born result, as a function of transverse momentum 
and pseudorapidity. Default parameters as specified in the text. 



p T ;> 2 GeV, where the scale and mass dependence is modest and the theoretical predictions 
appears to be well under control. The reduced theoretical uncertainty is a consequence both 
of the increase of the hard scale and of the increasing importance of the direct contribution, 
which receives only small radiative corrections above p T <; 2 GeV. This fact, together with 
the behaviour of the if-factor previously discussed, leads us to the conclusion that charm 
production in photon-photon collisions is under better theoretical control than in photon- 
hadron or hadron-hadron collisions. 

The dependence upon the parton density sets is mild and visible in the low-p T region 
only, where the resolved contribution is still sizeable (this is also illustrated in table |l|). The 
average Bjorken-x probed is of the order of 2m T /y&^, so that the photonic parton densities 
can in principle be tested in a region where they are poorly constrained by available data. 
Unfortunately, the sizeable theoretical uncertainty due to other sources in the cross section 
prediction at small p T prevents a discrimination of the different PDF sets. 

Finally, in the last plot of fig. ||] we show the sensitivity of the p T spectrum to the shape 
of the fragmentation function by varying the e parameter that enters the Peterson et al. 
form. We have considered, together with our default value e = 0.035, the more extreme 
choice e = 0.02, as obtained in recent fits [14] where the resummation of large collinear 
logarithms is taken into account. We remind the reader that the e values obtained from 



older fits [25| were significantly higher (e ~ 0.06), resulting in a softer p T spectrum. There 
seems to be a clear indication now that the degradation of the momentum of the parent 
charm quark should be smaller (the very same effect is also emerging in b physics). We also 
display the purely perturbative result without fragmentation, corresponding to e = 0. 
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da (eV -> e + e" D*X) / dp T [pb/GeV] 

Vs= 189GeV; |Tl(D*)|< 1.5 
m c =1.5 GeV; li R =\l t /2=m T ; 6=0.035 



GRV HO 
AFGHO 
GRS HO 




upper curves: NLO 
lower curves: LO/5 
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da (eV -> e + e" D*X) / dp T [pb/GeV] 
Vs= 189 GeV; fr|(D*)| < 1.5 
m c =1.5 GeV; 0^=11/2=11^; GRS 

e = (no FF) 
e = 0.02 
e = 0.035 




Figure 4: Dependence of the p T spectrum of the D* upon renormalization scale, 
charm mass, photonic parton densities and the parameter e of the Peterson et al. 
fragmentation function. Both the NLO (upper curves) and the LO (reseated by 1/5; 
lower curves) results are presented. 



From the figure one can conclude that the p T spectrum of the D* in 77 collisions is not 
sensitive enough to the value of e to discriminate between e = 0.035 and e = 0.02. On the 
other hand, the purely perturbative prediction without fragmentation function results in 
a much harder spectrum, which, as will be shown later in this section, is not favoured by 
the experimental data. Let us mention that the need to include a fragmentation function 
is not always obvious from a comparison with data; measurements of B meson production 
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do (eV -> eV D*X) / dp T [pb/GeV] 

Vs=189 GeV;|Ti(D*)|<1.5 

Hj^m,,; e=0.035; f(c->D*)=0.233; GRS 

m t =l .2 GeV; |j, R =2m T (dir) / tHj/2(les) 

m =1,5 GeV; fX R — it^ 

m c =1.8 GeV; fl^nL/X&T) I 2m T (res) 




- data: OPAL 
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da (eV ^ eV D*X) / di\ [pb] 

Vs = 189 GeV; 2 GeV < p T (D*) < 15 GeV 
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m c =1.2 GeV; p, R =2m T (dir) / metres) 

m =1.5 GeV; H- R =m T 

m c =1.8 GeV; (i^m^fdir) / 2m T (res) 
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Figure 5: Comparison between the NLO theoretical prediction and the OPAL data. 
The theoretical band is obtained by varying both the renormalization scale and the 
charm mass. The theoretical curves include the probability for a charm quark to 
fragment into a D* meson f(c — » D* + ) = 0.233 and take into account the fact that 
OPAL counts both the heavy quark and the heavy antiquark in their cross-section 
determination. 



in pp collisions at the Tevatron are consistent with a purely perturbative prediction, thus 
(provocatively) suggesting e = 0. 

The detailed studies of charm production at fixed-target experiments imply the relevance 
of non-perturbative effects, such as the intrinsic transverse momentum of the incoming 
partons fl]]. We attempted to model such effects by supplying the initial-state partons 
with a Gaussian-distributed transverse momentum. Although such a procedure is not fully 
consistent when radiative corrections are included, it can nevertheless provide an estimate 
of the relevance of these non-perturbative phenomena. We have studied this issue also in 
the case of photon-photon collisions and find that intrinsic transverse momentum effects 
are completely negligible, even if the average transverse momentum of the partons is chosen 
as large as 2 to 3 GeV. 

In fig. [5] we finally compare our predictions for the D* transverse momentum and pseu- 
dorapidity distribution with the OPAL measurement. Three different theoretical curves are 
shown, where the charm quark mass and the renormalization scale are varied as indicated 
in the figure, while the photonic parton density and charm-to- D* fragmentation function 
are kept fixed at the default choice. The theoretical results include the probability for a 
charm quark to fragment into a D* meson f(c — > D* + ) = 0.233 and take into account the 
fact that OPAL counts both the heavy quark and the heavy antiquark in their cross-section 
determination. We can conclude that the data agree reasonably well with the NLO pre- 
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dictions. More specifically, we can see that the shape of the distributions is reproduced 
quite well, while there is a small discrepancy in absolute normalization, of the same size 
as that observed when discussing the total rates. Of course, a definite conclusion will only 
be possible after the statistical errors affecting the data will decrease; at present, the data 
seem to suggest that the input parameters should be taken such as to yield the largest 
possible cross section. However, it is difficult to disentangle the various effects that could 
play a role here. We could indeed enlarge the cross section by taking a small mass value, 
a small renormalization scale, photonic densities with much softer gluon, or a combination 
of these three effects. Furthermore, the inclusion of NNLO is expected to further increase 
our prediction. In general, we can observe that theory undershoots the data when central 
values for the input parameters are adopted, as is the case for many such comparisons for 
charm production at various colliders. 

4. SUMMARY 

We have performed a full NLO QCD study of total and differential D* production rates 
in two-photon collisions at LEP2. Compared with charmed-meson production in hadron- 
hadron or photon-hadron collisions, the two-photon cross section appears to be under better 
theoretical control, mainly because of the dominance of the direct channel, where QCD 
corrections are smaller than in the case of the resolved channels. This consideration holds 
in particular for the moderate- and high-p T region and for the central rj region, which are 
directly probed by LEP experiments. We find the largest theoretical uncertainties, due to 
charm mass and renormalization scale dependence, at small p T 's, similarly to what has been 
observed in photon-hadron and hadron-hadron collisions. In this region, the contribution 
of the resolved channels is sizeable, and the theoretical predictions are mildly dependent 
upon the choice of the photonic parton densities. However, this dependence is smaller than 
that due to mass or renormalization scale choice, and this prevents a discrimination of the 
different density sets, even in the ideal case in which the charm mass would be known 
exactly. As far as the comparison with the data available at present is concerned, we find a 
good agreement in the shape of the p T and r/ distributions, and a reasonable agreement as 
far as the absolute normalization is concerned. In general, the comparison between NLO 
QCD and the experimental results reproduces the pattern already known from other types 
of collisions; namely, theory somewhat underestimates the data, and a special tuning of 
the input parameters is needed to improve the agreement. The mild dependence of the 
theoretical predictions for the visible cross section upon the input parameters appears to 
imply that it will be difficult to constrain any of these parameters by comparing theory 
with data. On the other hand, when the statistical significance of the measurements will be 
increased, charm production in photon-photon collisions will be a valuable tool in testing 
the underlying production dynamics. 
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APPENDIX: PHOTON-PHOTON CROSS SECTIONS 



In this appendix we describe some technical aspects of the computation of charm cross 
sections in photon-photon collisions. In particular, we will discuss in some detail the issue 
of scheme dependence of the direct and resolved production mechanisms. 

We consider the process 

e + + e~ — ► e + + e~ + 7 + 7 

Uq + q. 

The final-state positron and electron are scattered almost collinearly to the beam line; 
thus, the photons eventually producing the heavy quark pair are (almost) on-shell. The 
cross section can therefore be written as follows 

da ee (P e + , P e -) = J dx e +dx e -f!f\x e +)f!f\x e -)da(x e +P e +,x e -P e -), (5) 

~ (e) 

where da is the cross section for the process 77 — * QQ, and / 7 is the Weizsacker-Williams 



function [15|. We introduce the short-hand notation 

da ee = f\ e Kda*f\ e \ (6) 

which defines the * symbol. As mentioned in section 2, three mechanisms contribute to the 
two-photon cross section: 

da = da 7 -y + da\ T + da 2r , (7) 

where the three terms in the RHS of this equation denote the direct, single-resolved, and 
double-resolved components, respectively. Using the factorization theorems in QCD, the 
single- and double-resolved cross sections can be written as the convolution of the partonic 
hard-scattering cross sections da^, da^i and d&ij with the parton densities in the photon 

f (7). 
J i 

da lr = / i (7) *d* i7 + dfr 7 i*// 7) , (8) 

da 2r = // 7) *^*/j 7) , (9) 

Here, a summation over repeated indices is understood, and the indices i,j run over all the 
parton flavours u,u,- ■ ■ ,g. We can summarize the content of eqs. (0), @ and (||) by writing 

da = F^Kda ab *F^\ (10) 



13 



Direct 


Sinj 


de-res. 


Double-res. 


2^2 77 -> QQ 


70 ■ 


-> QQ 


qq — > QQ 








QQ 


2^ 3 77^QQff 


79 " 


-* QQ# 


55 <9<9s 




79 - 


^QQq 


qq -> QQ5 








95 <2<2<? 



Table 2: Partonic subprocesses contributing to heavy-quark production in 77 ces- 
sions next-to-leading order. 



where the indices a and b run over all parton flavours, plus the photon. The generalized 
densities in the photon are given by 

FW(x) = S aj S(l - *) + f^W (1 " <W) • (11) 

The three quantities on the RHS of eq. ([?]) are separately divergent order by order in 
perturbation theory, and scheme-dependent. The only quantity that is finite and scheme- 
independent at all orders is da. However, as mentioned in the introduction, the reaction we 
study has the special feature that the scheme dependence of the direct component and of 
the sum of the resolved components only starts beyond NLO. Although that still does not 
make these cross sections physically observable, it does give us a fairly clear idea of how 
often the initial photon fluctuates to a hadronic state at the energies considered. Let us 
therefore clarify the issue of scheme dependence by working out the specifics in some detail. 
We write the perturbative expansions of the partonic cross sections as: 

00 

d(Tryy ^eill ^ > Otgd(J^^j , (12) 

n=0 

oo 

dai-y = a eUi a s ^2a n s da^\ (13) 

n=0 

oo 

da l3 = a 2 s Y,* n sda%\ (14) 

n=0 

The expansion of da^i is identical to that of da^ in eq. (|13|). In what follows, we restrict 
ourselves to the first non-trivial QCD corrections (which amounts to retaining the first two 
terms in the sums of eqs. (0)^(0))' this * s phenomenologically relevant, since no complete 
NNLO result is available for heavy-flavour production. At this order, only (tree-level and 
loop) 2^2 and (tree-level) 2 — > 3 partonic processes contribute to the result; they are 
listed in table |2] for the three production channels. 

A general two-photon cross section such as da, being a physical observable, must be 
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scheme-independent. Therefore, an expression equivalent to eq. (|i"0|) is 

da = F^' *da' ab *F^\ 



(15) 



where the prime denotes that the generalized densities and the partonic cross sections, which 
are non-physical, are calculated in a scheme different from that used in eq. (|l0|). In QCD, 
the parton densities expressed in two different schemes can always be related as follows: 



F^' = 5 ai 5 + /W' (1 _ 5 ai ) , f™' = ft <] + a em Hi + a s K VJ ® &<\ 



(16) 



where Hi and are functions and distributions, respectively, with an expansion in a s 
starting from oP s . The symbol (g> denotes the convolution integral 



^ (g> B{z) = B (g> .A(z) = y dxdy <5(z - xy) A(x)5(z). 

Substituting eq. (jltf ) into eq. ([l5|), and using eqs. (12)-(|li|) we find: 

d<r = a 2 cra da^ + «em a s + Hi* da£j + d<r^ * Hi 



(17) 



~r Cl em Q!< 



(7)' 



^*d^'+d^' 



(7) 



+ * daf * + <*» f™ * dfi". 



.7 



p(7) 



.(1)' 



(7) 



+ (K ik 4 7) ) * d^f * ^ 7) ] + «emc4, « 



(18) 



Using only the first term in the expansion in a s of -ffj and Ku (which we denote as H\ 
and Kj-j\ respectively), and comparing eq. (|i~8| ) with eq. (|7|), we get 



(o) 









= < o) 


,.(0)' ,.(0) 

» d(T \i =d(?\i 




= <4V 


-nr. 






< r 






< - 












n k * a(J ki ! 








<- 


K (0) * d<T (0) 



da. 



(0) 



where use has been made of the fact that 

(A® B) * C = A* (B * C) 



(19) 
(20) 
(21) 
(22) 
(23) 

(24) 



Equation. (19) is nothing but the formal statement that, at leading order, the cross sections 



are independent of the scheme adopted. This equation also implies that, at this order, the 
direct, single-resolved, and double-resolved terms are separately physically meaningful. The 
inclusion of radiative corrections, however, changes the situation, as can be clearly seen from 
eqs. (|20|)-(p3|). For instance, when changing the scheme, the single-resolved cross section 
(eqs. (|H) and (|D) receives a contribution from the double-resolved part. In principle, this 
is also true for the direct cross section (eq. (ptj))- Note, however, that for heavy flavour 
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production, if i is a gluon, then H\ is zero, while if i is a quark, then da\^ is zero (this 
is not true for - say - jet production, where daffl is non- vanishing) . Therefore, at next-to- 
leading order the direct term is scheme- independent, while the single- and double-resolved 
terms are closely related, since one is entitled to add to them a finite piece without changing 
their sum. Clearly, when going to next-to-next-to-leading order, the direct cross section will 
also become scheme-dependent. 
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